Spin-phonon coupling in single Mn doped CdTe quantum dot 
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The spin dynamics of a single Mn atom in a laser driven CdTe quantum dot is addressed theoreti- 
cally. Recent experimental results^-^show that it is possible to induce Mn spin polarization by means 
of circularly polarized optical pumping. Pumping is made possible by the faster Mn spin relaxation 
in the presence of the exciton. Here we discuss different Mn spin relaxation mechanisms. First, 
Mn-phonon coupling, which is enhanced in the presence of the exciton. Second, phonon-induced 
hole spin relaxation combined with carrier-Mn spin flip coupling and photon emission results in Mn 
spin relaxation. We model the Mn spin dynamics under the influence of a pumping laser that injects 
excitons into the dot, taking into account exciton-Mn exchange and phonon induced spin relaxation 
of both Mn and holes. Our simulations account for the optically induced Mn spin pumping. 
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I. INTRODUCTION 

The tremendous progress in the miniaturization of 
electronic devices has reached the point that makes it 
crucial to address the effect of a single dopant in a de- 
vice and motivates the study of a single dopant spin to 
store digital information^. The manipulation of a single 
atom spin in a solid state environment has been demon- 
strated using several approaches, like scanning tunneling 
microscope on magnetic adatoms 5 ' 6 , or optical probing 
of NV centers in diamond^ and single magnetic atoms 
in semiconductor quantum dots, the topic of this paper. 
Single quantum dots doped with a single Mn atom can be 
probed by means of single exciton spectroscopy in pho- 
toluminescence (PL) experiments. This has been done 
both in H-Vlln2«fir-1£ and III-\2ML materials. In the 
case of single Mn doped CdTe dots, information about 
the quantum spin state of a single Mn atom is extracted 
from the single exciton quantum dot photoluminescence 
due to the one on one relation between photon energy 
and polarization and the electronic spin state of the Mn 
atomi£~— . This has made it possible to measure the 
spin relaxation time of a single Mn atom in a quantum 
dot under optical excitation, using photon autocorrela- 
tion measurements^, and to realize the optical initial- 
ization and readout of the spin of the Mn atomic—. 

The observation of Mn spin orientation under quasi- 
resonant optical pumping^— can be accounted for if the 
Mn spin relaxation time is shorter in the presence of a 
quantum dot exciton^^—. In that situation, resonant 
excitation of an optical transition associated to a given 
Mn spin projection results in the depletion of the laser 
driven Mn spin state, via Mn spin relaxation in the pres- 
ence of the exciton. Whereas theoretical understanding 
of the exchange couplings between electrons, holes and 
Mn spin in quantum dots permits to account for the ob- 
served PL spectr a 8 ' 10 ' 11 ' 18 , a complete understanding of 
the spin dynamics under the combined action of laser 
pumping, incoherent spin relaxation and coherent spin- 
flips is still missing. In this paper we make progress along 



this direction on two counts. First, we discuss different 
Mn spin relaxation mechanism, taking fully into account 
the interplay between incoherent dynamics due to the 
coupling to a reservoir and the coherent spin flips asso- 
ciated to exciton-Mn exchange in the quantum dot. Our 
calculations show that the most efficient Mn spin relax- 
ation channel, in the presence of the exciton, arises from 
a combination of phonon-induced hole spin relaxation, 
which turns the bright exciton into a dark, followed by re- 
combination to the ground enabled by dark-bright mixing 
due to Mn-carrier spin flip exchange. Thus, we provide a 
quantitative basis to a recently proposed scenario^. Sec- 
ond, we model the Mn spin dynamics with a rate equation 
for the Mn spin and the Mn plus exciton spin states that 
includes the spin relaxation rates between the few-body 
states calculated from microscopic theory. 

Our theory permits to model the experimental observa- 
tions and, importantly, it identifies the light-hole heavy- 
hole mixing as a crucial parameter that determines not 
only the PL lineshap o 10 ' 11 ' 18 but also the amplitude of 
several spin-relaxation mechanisms at play in this sys- 
tem. 

The rest of this manuscript is organized as follows. 
In section |TT] we review the Hamiltonian for a single 
Mn spin interacting with a single exciton in a quantum 
dot. The anisotropic Mn-hole coupling is derived from 
a simplifie d 18 ' 30 single-particle description of the lowest 
energy quantum dot hole states which affords analyti- 
cal expressions for the critical parameter in the theory, 
the light- hole heavy hole (LH-HH) mixing 30 . The de- 
pendence of the properties of the Mn-exciton states on 
the LH-HH mixing are discussed. In section (IIIII) we dis- 
cuss the Mn spin relaxation due to Mn-phonon coupling, 
both with and without an exciton in the quantum dot. 
Whereas this mechanism is probably dominant for the 
Mn spin relaxation in the optical ground state, it is not 
sufficient to account for the rapid Mn spin relaxation in 
the presence of the exciton. This leads us to consider 
other spin relaxation mechanisms. In section (|IV[) we de- 
scribe the spin relaxation of holes due to their coupling to 
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acoustic phonons, using a Bir-Pikus Hamiltonian. Using 
the simplified description of hole states, we obtain ana- 
lytical results for the hole spin lifetime in a non-magnetic 
dot, which are in agreement with previous work using a 
more sophisticated description of single hole states 3 -!. We 
then compute the lifetime of the exciton-Mn states due 
to hole spin relaxation. In section (|VI[) we present our 
simulations of the optical pumping process, using rate 
equations for the exciton-Mn quantum states, including 
the laser pumping, the spontaneous photon emission, the 
Mn and hole spin relaxation due to phonons. Our simu- 
lations account for the optical initialization and readout 
observed experimentally. 

II. EXCITON-Mn HAMILTONIAN 

In this section we describe a minimal Hamiltonian 
model that can accounts for the PL spectra of single Mn 
doped CdTe quantum dots. For that matter we need to 
consider both the Mn spin in the unexcited crystal and 
the Mn spin interacting with a quantum dot exciton. The 
peaks in the PL spectra are associated to the energy dif- 
ferences between the states of the dot with and without 
the exciton. 



A. Mn spin Hamiltonian 

Mn is a substitutional impurity in the Cd site of CdTe. 
Thus, it has an oxidation state of 2+, so that the 5 d elec- 
trons have spin S = |, resulting in a sextuplet 33 whose 
degeneracy is lifted by the interplay of spin orbit and the 
crystal field. In an unstrained CdTe, the crystal field 
has cubic symmetry which should result in a magnetic 
anisotropy Hamiltonian without quadratic terms. Elec- 
tron paramagnetic resonance (EPR) in CdTe strained 
epilayers^ show that Mn has an uniaxial term in the spin 
Hamiltonian. In a quantum dot there could be some in- 
plane anisotropy as well, which lead us to consider the 
following Hamiltonian: 

Ho = DMl + E{Ml - M 2 ) + gl i B B ■ M (1) 

where M a are the S — | spin operators of the electronic 
spin of the Mn. The eigenstates of this Hamiltonian are 
denoted by 4> m 

H \(f> m ) = E m \<f> m ) = E m £ <j> m {M z )\M z ) (2) 

where \M Z ) are the eigenstates of M z . In this paper 
we neglect the hyperfine coupling to the / = | nu- 
clear spin, which could affect the decay of the electronic 
spin coherence 1 -. The magnetic anisotropy parameters, 
E and D can not be inferred from PL experiments, 
which are only sensitive to the Mn-exciton coupling. 
EPR experiments 34 in strained layer could be fit with 
D = 12/zeV, E = and g = 2.0. Thus, the ground state 



should have M z = ±5, split from the first excited state 
by 2D. At 4 Kelvin and zero magnetic field, we expect 
all the six spin levels to be almost equally populated. We 
refer to these six states as the ground state manifold , in 
contrast to the excited state manifold, which we describe 
with 24 states corresponding to 4 possible quantum dot 
exciton states and the 6 Mn spin states. 

B. Single particle states of the quantum dot 

We describe the confined states of the quantum dot 
with a simple effective mass model. In the case of the 
conduction band electrons, we neglect spin orbit coupling 
and we only consider the lowest energy orbital, with wave 
function ^0(^)1 which can accommodate I electron with 
two spin orientations. 

In the case of holes, spin orbit coupling lifts the sixfold 
degeneracy of the top of the valence band into a J = § 
quartet and a J = \ doublet which, in CdTe, is more than 
0.8eV below in energy. Confinement and strain lift the 
fourfold degeneracy of the J = | hole states, giving rise 
to a mostly J z — ±| heavy hole doublet and a almost 
J z = ±i light hole doublet. Importantly, it is crucial 
to include LH-HH mixing to describe the experimental 
observation. 



1. Effect of confinement 

The top of the valence band states are described in 
the kp approximation with the so called Kohn-Luttinger 
Hamiltonia n 35 ' 36 . For that matter, the crystal Hamilto- 
nian is represented in the basis of the 4 topmost J = | 
valence states of the T point. We label them by J = 
I , J z . The resulting kp Hamiltonian can be written as 

Wholes = %KL 

Hkl- V^JiJjhkj+KHsJzB (3) 

where ki are the wave vectors, Ji are the spin | matrices, 
and V^j h are coefficients given in the appendix IA1I The 
last term accounts for the Zeeman coupling to an external 
field along the growth direction. 

The kp Hamiltonian for states in the presence of a 
quantum dot confinement potential that breaks transla- 
tional invariance reads: 

n kp = -h 2 J2 y^JiJA d 3 + wbJzB + v(r)s JzJ , 

(4) 

In general, the numerical solution of equation can 
be very complicated. Following previous wor k 18 ' 30 , we 
make two approximations that permit to obtain analyt- 
ical solutions. First, we model the quantum dot with 
a hard-wall box-shape potential. The dimensions of the 
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box are L x , L y and L z . This permits to write the wave 
function as a linear combination of |J — - 
multiplied by the confined waves 



^) d z) states 
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(5) 

Our second approximation is to restrict the basis set to 
the fundamental mode only, = 1. As a 

result, the quantum dot Hamiltonian reads 



H 



k P 



where 



(didj) = J ^1,1,1 (f)di9j-Vi,i,i(r) = % f^? 



(0) 



(7) 



Thus, within this approximation, the quantum dot hole 
states are described by a 4*4 Kohn Luttinger Hamilto- 
nian where the terms linear in fc.; vanish and the kf terms 

are replaced by (jf^j ■ The resulting matrix Hconf has 

two decoupled sectors denoted by + (J z = +§ ,J Z = —5) 
and — (J z = 
basis we have: 
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where b = kvbB and P, Q and i? are given in the ap- 
pendix. In order to find the corresponding energies and 
wavefunctions it is convenient to write these matrices as: 
%± = a± + h± ■ a where a are the Pauli matrices act- 
ing on the pseudospin | space of the + and — spaces, 
a± = P =F b/2 and 

h± = (R,0,QTb) = \h±\ (Sin0±,O,Cos9±) (11) 

We keep only the two ground states (heavy- hole like), 
denoted by | ff) and | JJ-)? which are several meV away 
from the light-hole like states. The ground state doublet 
for the quantum dot holes states so obtained, neglecting 
strain, can be written as 



ft) 
ft) 



2)+sm-\--) 

-3. . 0-. 1. 
->+- T l + 2> 



fl2) 



Thus, the LH-HH mixing parameters, depend on the 
dot dimension, L iy on the Kohn Luttinger parameters, "fi 
and on the applied magnetic field B. 



2. Effect of homogeneous strain 

We now consider the effect of the strain that arises from 
the lattice mismatch between the CdTe quantum dot and 



the ZnTe substrate on the J 



states of the valence 



band. It has a similar effect that confinement, resulting in 
a splitting of the J — | manifold and a mixing of the LH 
and HH states. The Hamiltonian that describes the effect 
of strain, as described by the strain tensor , on the top 
of the valence band states in zinc-blende semiconductors 
was proposed by Bir and Pikus. We can write the Bir 
and Pikus (BP) Hamiltonian as^ : 

( % \ / 
rtBP = I a - -r J (e xx + e yy + e zz ) + 

+36 J? en + V3d((J x J y + J V J X ) e xy + c.p.) (13) 



where c.p. stands for cyclic permutation, and a — 
-3.4eV, b = -1.2eV, d = -5.4eV for CdTe^. 

For CdTe quantum dots grown in ZnTe, we mainly 
consider the effects of strain anisotropy in the growth 
planei 3 - and describe the strain by the average values 



of 



xy 



and 



In this approximation the BP 



Hamiltonian is reduced to a block diagonal matrix in the 
(+§,-5, +5,-3) basis: 



H 



BP 
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(14) 



(15) 



(16) 



where A;^ is the HH-LH splitting, p s the strained in- 
duced amplitude of the HH-LH mixing and <p s the angle 
between the strained induced anisotropy axis in the quan- 
tum dot plane and x (100) axis. Importantly, the effect 
of confinement and the effect of strain have a very simi- 
lar mathematical structure. They both split the LH and 
HH levels and mix them. The main difference lies in the 
mixing term, which is real for the confinement Hamil- 
tonian controlled by the shape of the quantum dot and 
complex for the BP Hamiltonian depending on the strain 
distribution in the quantum dot plane. 



3. Combined effect of confinement and strain 

We finally consider the combined action of confinement 
and strain described by H holes = Hccmf + HBP- Summing 
the Hamiltonians of equations ((5J and (p~4|) to obtain two 
decoupled matrices for the + and — subspaces. They can 
be written as 



n to t,± = A ± +H±-a 



(17) 
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where A± = P =p 



and 



R + p s cos(2ip s ), ±p s sin(2ip s ), QTb ^- ) (18) 



It is convenient to express the ground state doublet asso- 
ciated to "Hhoios hi terms of the spherical coordinates of 
the vectors H±, \H±\, 0± and 



Cos- 



Cos- 



-1, 



where 
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R + p s e ±2ilfi ° 
\R + p s e ±2i f>° 



(19) 
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Expectedly, this expression is formally very similar to 
that of equation (|12l) . 

Formally, we express eq. (fT9|) as 



(21) 



Wh) =^2c h (j z )\j z 



dependent constant j3 and on a quantum dot dependent 
property, the probability of finding the hole at the Mn 
location. 

After a straightforward calculation we obtain the effec- 
tive Mn-hole coupling spin model working in the space 
(M,ah) of dimension 12 as a function of the hole wave 
function parameter 9: 



V,, — Mn — Jhx 

M X <7 X + jhyMyCTy + jhzM z a 



(24) 



where the jhi are dimensionless coefficients given, for B 
and 9 = 9+ = 6L by: 



Jhx 
3hy 
Jhz 



;]h_ 
2 



U/ZSine + 1 - Cos9 
y {yiSin9 - 1 + Cos9 
y(l + 2Cos9) 



(25) 
(26) 
(27) 



Notice that for 9 = there is no LH-HH mixing and 
we have jhx — jhy = and jhz = \jh- In this extreme 
case the Mn-hole coupling is Ising like and M z and a z are 
conserved. This limit is a good starting point to model 
hole-Mn coupling in CdTe quantum dot a 18 ' 39 



Both in equations (|T2|) and (fT9l) the (+§,—5) sector 
is decoupled from the (— 1>+|)- Whereas, this is not 
true in general, it is sufficient to account for the correct 
symmetry of a variety of exchange couplings between the 
hole and both the Mn and the electrons. 



C. Effective Mn-carrier exchange Hamiltonian 

1. Hole-Mn Hamiltonian 

We now consider the exchange coupling of hole spin (J) 
and the Mn spin (M). The leading term in the exchange 
interaction is the Heisenberg operato r 33 ' 38 : 



Voxch = I38(f h - r M )J ■ M 



(22) 



where (3 is the hole-Mn exchange coupling constant. For 
Mn in CdTe we have /3N =0.88 eV, where N is the 
volume of the CdTe unit cell^ 3 .. The exchange interaction 
is taken as short ranged, the Mn atom is located at tmu 
and J are the spin | angular momentum matrices. We 



represent the operator (j2"2")l in the product basis \M) x a^. 
Thus, the exchange operator in the product basis reads: 



(M\(ah\V cxch \M')\a' h ) = 
P\1>o (r Mn )\ 2 J2 {M\M a \M') I (a h I JaW h ) 



(23) 



2. Electron-Mn Hamiltonian 

In analogy to the hole-Mn bare coupling, the electron- 
Mn coupling reads: 



V e - 



Mu 



aS{r e - r M )S ■ M 



(28) 



where S is the spin of the electron. Since the spin orbit 
coupling has a very small effect on the s like conduction 
band, the effective exchange for the quantum dot electron 
and the Mn is also a Heisenberg term given by 



V e - 



: Je S ■ M = J. (S X M Z + SyMy 



S Z M Z ) (29) 



where j e = a\ipo(^Mn)\ ^ s electron-Mn coupling 
constant, which depends both on a material dependent 
constant a and on a quantum dot dependent property, 
the probability of finding the electron at the Mn loca- 
tion. In our model we take the same orbital wave func- 
tion for the confined electron and hole, so that the ratio 
je./ jh — ct/P — 4 for CdTe. A deviation from this ex- 
pected ratio is indeed observed in real CdTe quantum 
dots£. 



D. Exciton-Mn wavefunctions and energy levels 

1. Hamiltonian 



where ij^oif) is the envelope part of the heavy hole wave 
function, eq. ([S]), and jh = 0\i>o(rMn)\ 2 is the hole- 
Mn coupling constant, which depends both on a material 



The effective Hamiltonian for the exciton in a single 
Mn doped CdTe quantum dot is the sum of the single 
ion magnetic anisotropy Hamiltonian, the Mn-electron 



and Mn-hole exchange coupling and the electron-hole ex- 
change coupling 



where 



T-L = Hs + Ve-Mn + Vh-Mn + Ve-h 



Ve-h — jehS z (J h 



(30) 



(31) 



is the electron hole exchange coupling, neglecting trans- 
verse components. Electron hole exchange is ferromag- 
netic (j e h < 0) and splits the 4 exciton levels into two 
doublets, the low energy dark doublet ("fit, -H-Di denoted 
by X — ±2 and the high energy bright doublet (fti, JJ-t) 
(X = ±1). 

Since we consider two electron states (S z =t,l), 
two hole states (ah =1)", JJ-), and six Mn states M z = 
±|,±|,±i, the Hilbert space for the Mn-exciton sys- 
tem has dimension 24. Whereas we do obtain the exact 
eigenstates of Hamiltonian (fBTJf by numerical diagonal- 
ization, it is convenient for the discussion to relate them 
to eigenstate of the Ising, or spin conserving part, of the 
Hamiltonian: 

T~L = Rising + %flip (32) 

where 

•Hi^ng = DMl + j eh S z a h + j e S z M z + j h M z a h (33) 

and 

H fli p = E{Ml - M^y) + j e (S X M X + SyMy) + 

+ (jhx cr x M x + j hy <Ty My ) (34) 

If we expand jh x and jh y in the series of LH-HH mixing 
parameter 8, they are the same in the first order of 9. 
For simplicity, we take 



]h± = Jhx = Jhy = Jh 



(35) 



in the following calculation. In the case of a LH-HH 
mixing induced by the anisotropy of the confinement de- 
scribed by a hard-wall box shape potential, we get from 
the Kohn-Luttinger Hamiltonian: 



(36) 



The eigenstates of "Hising are trivially given by the 
product basis 

\P) = \M z )\S z )\a h ) (37) 

with eigenenergies: 

E P = DM? + iehS z <r h + j e S z M z + Jh M z a h (38) 




FIG. 1: (Color online) Scheme of the energy levels of the 
quantum dot exciton interacting with 1 Mn when spin-flip 
terms are neglected. 



Since the magnetic anisotropy term DM^ is present both 
in the ground state and exciton state manifolds, it does 
not affect the PL spectra of the bright excitons. Within 
this picture, for each of the 6 possible values of M z , there 
are 4 exciton states. We use a short-hand notation to 
refer to the Ising states Px (M z ) where X = ±1, ±2 labels 
the spin of the exciton, X — S z +a z . An energy diagram 
for the exciton levels, within the Ising approximation, is 
shown in figure |T]). 

The PL spectra of a single Mn doped quantum dot 
predicted by the model of Ising excitons, ie, neglecting 
the spin flip transitions, features 6 peaks corresponding 
to transitions conserving M z . For the recombination of 



er+ excitons (S z = 
corresponds to M z 

M z 



j , ct/j =ff) the high energy peak 
+ | and the low energy peak to 
■| on account of the antiferromagnetic coupling 
between the hole and the Mn. In the case of a~ excitons 
the roles are reversed, but the PL spectrum is identical 
at zero magnetic field. 



2. Wave functions 

When spin-flip terms are restored in the Hamiltonian, 
the P states are no longer eigenstates, but they form a 
very convenient basis to expand the actual eigenstates of 
H, denoted by |\P„): 



|* n )=£* n (P)|P) 
p 



= V n {X,M z )\X,M z ) (39) 



X,Mz 



In most cases, there is a strong overlap between ty n and 
a single state \P). This is expected for several reasons. 
First, the single ion in plane anisotropy is probably much 
smaller than the uniaxial anisotropy, D E. Second, 
the electron-hole exchange, which is the exchange energy 
in the system, splits the dark and bright levels. Thus, 



6 




0.04 



FIG. 2: (Color online) Left panel: Evolution of the exciton 
levels as a function of the LH-HH mixing parameter Jh±- 
Right panel: Evolution of the IPZ as a function of the LH- 
HH mixing parameter. The inset presents the evolution of 
the energy for all the 24 exciton levels. 



both electron and hole spin flip due to the exchange with 
the Mn spin is inhibited because they involve coupling 
between energy split bright and dark excitons. In addi- 
tion, the electron Mn exchange is smaller than the hole 
Mn exchange, whose spin-flip part is proportional to the 
LH-HH mixing and approximately 10 times smaller than 
the Ising part. In order to quantify the degree of spin 
mixing of an exact exciton state ty n , we define the in- 
verse participation ratio: 



IPZ„ 



p 



\^n(P)\ 4 



(40) 



This quantity gives a measure of the derealization of 
the state \I/ n on the space of product states of eq. (1571) . 
In the absence of mixing of different P states, we have 
IPZ n = 1. In the case of a state equally delocalized in 
the 24 states of the P space, we would have $? n (P) = — j= 

and IPZ n = Jj. 

In figure © we plot the evolution of both the energy 
(left panel) and the IPZ (right panel) as a function of 
JhX> the LH-HH mixing parameter, of four states de- 
noted by their dominant component at Jh± = 0. For our 
choice of exchange constants, two of them | + 1, — |) and 
| — 2, — A) are almost degenerate at J^± = 0, which means 
that the hole-Mn exchange compensates the dark-bright 
splitting, and couple these states via a hole-Mn spin flip. 
As a result, their energy levels split linearly as a function 
of Jh± and the wave-functions have a large weight on the 
two product states for finite Jhj_- In contrast, the other 
two levels shown in figure ©, | + 1, — |) and | — 2, — ~) 
are not coupled via hole-Mn spin flip. As a result, their 
energies shift as a function of J^j_ due to coupling to 
other states, and their IPZ undergoes a minor change, 
reflecting moderate mixing. 



3. Exchange induced dark-bright mixing 

The most conspicuous experimentally observable con- 
sequence of the exchange induced mixing, is the transfer 
of optical weight from the bright to the dark exciton, 
which results in the observation of more than 6 peaks in 
the PL. This can be understood as follows. The spin-flip 
part of the hole-Mn interaction couples the bright exci- 
ton | + 1, M z ) to the dark exciton | - 2, M z + 1). Thus, 
a state with dominantly dark character | — 2, M z + 1) 
and energy given, to first order, by that of the dark ex- 
citon, has a small but finite probability of emitting a 
photon through its bright component, via a Mn-hole co- 
herent spin-flip. Thus, PL is seen at transition energy 
of the dark exciton. Reversely, nominally bright excitons 
loose optical weight due to their coupling to the dark sec- 
tor. Importantly, the emission of a photon from a dark 
exciton with dominant Mn spin component M 2 , entails 
carrier-Mn spin exchange, so that the ground state has 
M z ± 1. 

According to previous theory worki£ the rate for the 
emission of a circularly polarized photon from the exciton 
state ^ n to the ground state 4> m reads: 



F = F 



>{M s )V> n {M X) X = ±\) 



where 



r 



47refi,c 3 



(41) 



(42) 



is the recombination rate of the bare exciton, uj is the 
frequency associated to the energy difference between the 
exciton state n and the ground state m, c is the speed 
of light, e is the dielectric constant of the material, d cv 
is the dipole matrix element. From the experiments, we 
infer Tq =0.5 ns _1 

In the absence of spin- flip terms, the matrix 
would have only non-zero elements for n = \X = ±1, M z ) 
states connected to m = M z states. The presence of 
spin-flip terms in the Hamiltonian enables the recombina- 
tion from exciton states with dominant dark component. 
In figure we represent the matrix elements r^ m /Fo 
for j e h = — 0.73meV, jh = O.'S&meV, j e = — 0.09meV, 
j h ± = 0.036meV", D = O.OlmeV and E = OmeV. It 
is apparent that the recombination rates from the dark 
states are, at least, 2 times smaller (a and b) than those 
of the bright states. For Tq= 0.4 ns _1 , the lifetime of 
the dark excitons (a and b) are in the range of 3 ns. 
Thus, this provides a quite efficient Mn spin relaxation 
mechanism, provided that a dark exciton is present in 
the quantum dot. 

The recombination rate matrix, together with the non- 
equilibrium occupation of the exciton states, P n , deter- 
mines the PL spectrum^ for ± circular polarization: 



E 



P Y ± 

r n>- n i 



,6 (tkj - E n - E m ) 



(43) 
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where 
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FIG. 3: (Color online) Recombination rates of the excitons 
levels in a Mn doped quantum dot hi/ho: a : | + 2, +| >— > 
\ + | >, b : | -2,-| >-*■ ] - | >, c : | - 1,+f | + | >, 
d : | + 1, — f >— > | — f >, states from e to n are bright 
excitons. Other states are mainly dark excitons with a small 
bright component. 



In a typical PL spectrum 3 , the dark peaks are, at most, 2 
times smaller than the bright peaks. Since T^ m is at least 
2 times smaller for dark states, this implies a larger oc- 
cupation of the dark states. Thus, we can infer a transfer 
from the optical ground state to the dark exciton states, 
via the bright exciton states. This transfer requires an 
incoherent spin flip of either the electron or hole. Below 
we show that phonon induced hole-spin relaxation pro- 
vides the most efficient channel for this bright to dark 
conversion. 



III. MN SPIN RELAXATION DUE TO 
SPIN-PHONON COUPLING 



In this section we discuss the Mn spin relaxation in the 
absence of excitons. In the absence of carriers and given 
the fact that Mn-Mn distance is comparable to the dot- 
dot distance (100 nm for a dot and Mn density of about 
10 10 cm -2 ), which makes direct super-exchange negligi- 
ble, the Mn-phonon coupling should be the dominant, 
albeit small, Mn spin relaxation mechanism. Transverse 
phonons induce local rotations of the lattice. Since the 
crystal field, together with spin orbit coupling, deter- 
mines the Mn magnetocrystalline anisotropy, the phonon 
induced lattice rotation^ acts as a stochastic torque on 
the Mn spin, resulting in spin relaxation. 

The atomic displacement at point r in the crystal is ex- 
pressed in terms of the phonon operators with wave vec- 
tor q, polarization mode A = T1,T2,L, frequency uj\{q) 
and polarization vector e\{q^: 



u(r) 



ff.X 



a iq-r 



(44) 



Ux{q) = 



2pu x {q)V 



(45) 



and V and p are the volume of the crystal and the mass 
density respectively. In a zinc-blende structure there are 
two transverse acoustic (TA) phonon branches and one 
longitudinal acoustic branch (LA) . Following Woods 3 -! we 
have: 



Stai = — (qxq z ,q v qz, -q±) 
qq± 



e T A2 = — {q y ,-qx,0) 
q± 



(46) 



(47) 



qy . These vectors satisfy 



where q = \q\ and q± = 

q- &TAi = 0, and erAi • &TAj = $ij The longitudinal mode 
has e L A = \q- 

The lattice rotation vector is given by^£ 

5$(r) = V x u{r) (48) 

so that only the transverse modes contribute. Within 
this picture, the Mn spin-phonon coupling can be written 



V M - P h = i n ,M -S^ifMn) 



(49) 



Without loss of generality we can set the Mn position as 
the origin, rMn = 0. Equation (14T)1) couples the Mn spin 
to a reservoir of phonons whose non-interacting Hamil- 
tonian is 



(50) 



q.X 



Within the standard system plus reservoir master 
equation approach, we have derived the scattering rate 
from a state n to a state n' , both eigenstates of the single 
Mn Hamiltonian Hq, due the emission of a phonon. In 
order to use a general result for that rate (|B7|) . derived 
in the appendix (jB]), we need to express the spin-phonon 
coupling (1491) using the same notation than in equation 
(EH): 



= i 2 Ux{q)f n ,n' ■ {q x e\(q)) 



where 



fn 



H ,M 



(51) 



(52) 



We compute now the scattering rate due to a single 
phonon emission assuming 3 dimensional phonons de- 
scribed above. The rate reads: 

r "--' = i^fc ( " s(A)+1) £ &.»(/£.»)• 

6,6'— a:, y,z 

(53) 

where c = 1.79kms is the CdTe speed of sound^, p = 
5870kgm~ 3 is the mass density of the CdTe unit celtf2 
and ng(A) = efsl l l _ 1 ■ The |A| 3 factor comes from the 
dependence of the phonon density of states on the energy. 



A. Mn spin relaxation in the optical ground state 

We now discuss the relaxation of the Mn electronic spin 
due to spin-phonon coupling without an exciton in the 
quantum dot. According to our experimental results*^ 3 -, 
the Mn spin relaxation time in our samples is at least 
5 /as. 

If we take E — 0, the transition rate between the ex- 
cited states \(j) n ) = \M Z = +§) and \cj) n >) = \M Z = + §), 
via a phonon emission, is given by: 



r , — 

± n—>rr 



640|L>| 5 
3nph 4 c 5 



(n B (A) + l) 



(54) 



The dependence on D 5 comes both from the density of 
states of phonons paw 3 and the square of the Mn 
phonon coupling, which is proportional to the anisotropy, 
and gives the additional D 2 factor. Whereas the uniaxial 
anisotropy of Mn in CdTe quantum wells has been de- 
termined by EPR^i, the actual value for Mn in quantum 
dots is not known and can not be measured directly from 
single exciton spectroscopy of neutral dots. Therefore, in 
figure (j4j we plot the lifetime for the transition of the Mn 
spin from | to | , due to a phonon emission, as a function 
of D. We take D in a range around the value reported 
for CdTe:Mn epilayers, D = 12^eV 34 . We find that the 
spin lifetime of Mn in the optical ground state can be 
very large. Even for D = 20fieV the Mn spin lifetime 
is in the range of 0.1 seconds, well above the lower limit 
for the Mn spin relaxation reported experimentally^ 3 -. 
Whereas we can not rule out completely that the Mn 
spin lifetimes that long, there are other spin relaxation 
mechanisms that might be more efficient that the Mn- 
phonon coupling considered above, like the coupling of 
the Mn electronic spin to nuclear spins of Mn and the 
host atoms 1 -. 

Since part of the A 5 scaling arises from the uj 3 scaling 
of the phonon density of states, we have explored the pos- 
sibility that phonons localized in the wetting layer could 
be more efficient in relaxing the Mn spin. For that mat- 
ter we have considered a toy model of two dimensional 
phonons confined in a slab of thickness L — 2 nm. The 
resulting Mn spin relaxation rate for those reads: 



r„ 



16h 3 c 4 pW 



b,b'— x,y,z 



(55) 

where W is the width of the sample and A is a diagonal 
matrix with A xx = 1, A yy = 1, A zz = 2. 

In figure Q we plot the associated spin lifetime in this 
case, taking W — 2nm and show how it is at least 100 
shorter than for 3D phonons, but still we would have 
Ti ~ 1 ms for D = 20/ieV. 



B. Mn spin relaxation in the presence of an exciton 

Here we discuss how the Mn spin relaxation due to 
Mn-phonon coupling is modified when an exciton is in- 
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FIG. 4: (Color online) Left panel: Lifetime of the (+| to 
+ |) transition in the optical ground state at zero field as a 
function of the magnetic anisotropy energy splitting D. Right 
panel: Lifetime of the same transition in the presence of a +1 
exciton for different values of D. The rates are calculated 
for a 3 dimensional (3D) and a 2 dimensional (2D) density of 
states of acoustic phonons. 



teracting with the Mn. The Mn-phonon coupling is still 
given by Hamiltonian (l49l) . with Ho given by eq. JT]). We 
assume that the only effect of the exciton on the Mn is 
to change the energy spectrum and mix the spin wave- 
functions, giving rise to larger spin relaxation rates, due 
to the larger exchange-induced energy splittings. 

In the presence of the exciton, the Mn-phonon coupling 
results in transitions between different cxciton-Mn spin 
states, n and n! . As we did in the case of the Mn without 
excitons, we need to express the spin-phonon coupling 
(|49| using the same notation than in equation (|Blj) . 

For that matter we define the matrix elements 



F n , n , = (*„i [h ,m\ i*; 

Y y n (X,M)*y n (X,M')*l 



M,M> 



X,Mz,M'. 



where 



f M ,M' = {M\ H ,M \M') 



(56) 



(57) 



M and M' stand for eigenstates of the Mn spin operator 
M z . Thus, in the exciton-Mn spin states basis, the Mn- 
phonon coupling reads: 



= i 2 Ux{q)F n ,n>{M) ■ (qx e x (q)) 



(58) 



Notice how if we neglect the spin mixing of the exciton 
states we have F n _ n > — /m,a/' and the only difference in 
the scattering rates arises from the larger energy split- 
tings in the presence of the exciton. 

Using the equation (|B7j) for the phonon induced spin 
relaxation rate, and in analogy with equation (|53[) we 



write: 



n—>n' — 



I2irph 4 c 5 



(n B (A) + l) J2 

b,b'= x,y,z 

(59) 

In figure ((4]) we see how Mn-phonon spin relaxation is 
much faster in the presence of the exciton. Ignoring the 
difference arising from the spin mixing, we can write the 
ratio of the rates as: 



Tru-n' (G) 



Ax 
A G 



(60) 



The energy splitting associated to the | to | spin flip 
in the ground state is AD. In the presence of the exci- 
ton the energy splitting of the same transition would be 
AD + jh — j e - If we take D = 12/^eV, jh = 360^eV and 
j e = — 90/ieV the ratio yields ~ 10 3 . From the experi- 
mental side we know that Tiq > 5/is and, in the presence 
of the exciton T\ ~ 50ns. Thus, the ratio could be ac- 
counted for by this mechanism. However, in order to 
have Tig = 5 /its we would need to assume an unrealisti- 
cally large value for D. Thus, we think that another spin 
relaxation mechanism must be operative in the system 
when the exciton is in the dot which makes it possible 
to control the spin of the Mn in a time scale of 50 ns. 
In the next sections we discuss the hole spin relaxation 
due to phonons as the mechanism that, combined with 
Mn-carrier exchange, yields a quick Mn spin relaxation 
in the presence of the exciton. 
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FIG. 5: (Color online) Left panel: Hole spin-flip rate as 
a function of the dot size and shape. L y is fixed at 6nm, 
scaning of L x change the LH-HH mixing. In the inset the 
ratio L y /L x is fixed at 1.2 so the shape of the dot are fixed, 
only the size of the dot changes. One can see that the hole 
spin-flip rate is a size sensitive quantum quantity, the rate is a 
semi-exponential function of the size of the dot. Right panel: 
Hole spin-flip rate as a function of the energy splitting for 
two different values of the quantum dot anisotropy, ie LH-HH 
mixing. 



effective hole-phonon Hamiltonian is obtained by project- 
ing the BP Hamiltonian onto this subspace: 



IV. HOLE SPIN RELAXATION IN NON 
MAGNETIC DOTS 



Vh-ph 



.n= £ f l3 h ^{q)\a h ){a' h \e %1 {q) (64) 

ij,fch,tr' h 



A. Hole-phonon coupling 

We now consider the relaxation of the hole spin due 
to hole-phonon coupling. We consider first the case of 
undoped quantum dots. The coupling of the spin of the 
hole to phonons can be understood extending the Bir 
Pikus Hamiltonian to the case of inhomogeneous strain 
associated to lattice vibrations: 



M 2 I fir, d n 



(61) 



(62) 



It is convenient to write the strain tensor field as: 

so that we write: 

(63) 

We consider the coupling of the ground state doublet, 
formed by states and 4J- 7 to the phonon reservoir—. The 



Here \<Jh) denotes the quantum dot state defined in eq. 
(fT9l) and the coupling constant reads 



T a h,<? h (-*. _ 



E 

3z ,3'z 



W = £ VijCWJCv&WMJjltilt (65) 



where Ig — J |-0o(r)| 2 e l? ' r c?r. Hamiltonian shows 
how the absorption or emission of a phonon can induce 
a transition between the two quantum dot hole states, ft 
and JJ-. 

We now calculate the time scale for the spin relaxation 
of a single hole in a non magnetic dot under the influ- 
ence of an applied magnetic field so that the hole ground 
state doublet is split in energy. In order to compute the 
transition rate for decay of the hole from the excited to 
the ground state we use again the general equation (|B7[) . 
For that matter, we express the hole-spin coupling (|6~4")) 
as: 



phonon 



q,Xa h ,o 



q,\ 



\<Th 



)K\(b\ q 



(66) 
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where 



h3 



(67) 



B. Calculation of hole spin- flip rates with simple 
model 



In order to illustrate the physics of the phonon-driven 
hole spin relaxation we consider the case of a single hole 
in a non-magnetic dot under the influence of an applied 
magnetic field. For that matter, we compute the Hamil- 
tonian (|6"T|) using the wave functions from the simple 
model of confined holes defined in eq. (|T2"1) . We focus 
on the non-diagonal terms in the hole spin index, i.e., 
the terms that result in scattering form ft to ft due to 
phonon emission. 

Importantly, the BP Hamiltonian couples hole states 
that differ in, at most, two units of J z . Thus, in the 
absence of LH-HH mixing, the BP Hamiltonian does not 
couple directly the ft and ft states. Transitions between 
ft and ft states, as defined in eq. (fT9)) . are only possible, 
through one phonon processes, through the e yz (J y J z + 
J z Jy) and e zx (J z J x + J X J Z ) terms in the Hamiltonian. 
After a straightforward calculation we obtain: 



(e y z(q) - ie zx (q)) (68) 



The important role played by the LH—HH mixing #1.2 is 
apparent. Using equation (|B7j) it is quite straightforward 
to compute the rate for the 3 phonon branches. They are 
all proportional to 



l 1r-^ _ lg7r u u' i)m 



A 3 

P h 4 c 5 



(69) 



with coefficients |, 



1 and § for the TA1, TA2 and L 



modes respectively. Here, D u i stands for the deforma- 
tion potential of Kleiner- Rothi 3 -, following reference 0, 
D u i — — 3 ^ d . p for the mass density of CdTe, c for its 
transverse speed of sound, and A for the energy split- 
ting between the ft and ft states, which is proportional to 
the external magnetic field B. In figure §5§ we plot the 
rates Ttai,^tA2^l, as well as their sum as function of 
the dot size (left panel) and as a function of the energy 
splitting between the initial and final hole state, A (right 
panel). We see how hole spin relaxation rates can be in 
the range of V ~ l/(40ns). 

The results of figure ((5|) suggest that for sufficiently 
high A, as those provided by the Mn-hole exchange, the 
hole spin can relax in a time scale of 30ns. These numbers 
are in the same range than those obtained by Woods et 
al 3 ^. As we discuss in the next section, these spin flips, 
together with Mn-carrier exchange, can also induce Mn 
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FIG. 6: (Color online) Scheme of the Mn spin flip chan- 
nels due to the combined action of hole-phonon coupling and 
carrier-Mn exchange. 



spin relaxation in a time scale much shorter than the one 
due to Mn-phonon coupling only. 

Importantly, the rate is finite only if B\ — 6*2 ^ 0, which 
is the case in the presence of an applied magnetic field. 
This indicates that, within the simple model of eq. (|12j) . 
the non-diagonal terms in the hole-phonon Hamiltonian 
vanishes identically. This is not a general feature 
of (f64|) . but rather an particular property of the simple 
model (fi"2"|) . In particular, the non-diagonal term in (fT2|) 
is non-zero at zero field as soon as the (+§,—5) states 
have also some weight on the +i components, which hap- 
pens both when a more realistic model for confinement 
is used or when homogeneous strain components e yz and 
e zx are included. 



V. SPIN RELAXATION IN MAGNETIC DOTS 
DUE TO HOLE-PHONON COUPLING 

The results of the previous sections indicate that, be- 
cause of their coupling to phonons, the hole spin life- 
time in a non-magnetic dot is much shorter than the Mn 
spin lifetime. Here we explore the consequences of this 
phonon-driven hole spin relaxation for the single exciton 
states in a dot doped with one magnetic atom. The lead- 
ing process results in a Mn spin conserving decay from 
the bright exciton to the dark exciton state, via hole- 
spin flip in a time scale in the 10 ns range. Combined 
with the optical recombination of the dark state, made 
possible via Mn-hole or Mn-electron spin flip, provide a 
pathway for exciton induced Mn spin relaxation in a time 
scale under 100 ns, as observed experimentally^—. 

We also explore the scattering between two bright 
states enabled by the combination of phonon induced 
hole spin relaxation and Mn-carrier exchange. The life- 
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times of these processes is in the range of 10 3 ns and 
higher, and therefore they are probably not determinant 
for the optical orientation of the Mn spin in the sub- 
microsecond scale. 



A. Exciton-phonon coupling in magnetic dots 

The Hamiltonian that couples the exciton states 
to the phonons is derived by projecting the hole-phonon 
coupling Hamiltonian (|64|) onto the exciton states ()39|) . 
The result reads: 

V X _ phon = £ \*n){*n>\V^' (b{ q + b^ q ) (70) 

n : n'q\ 

where 

V?,a' = E V^<V n (M z ,X)K'(M z ,X') (71) 

M z ,<7e,a h ,a' h 

where X — {a e ,<Jh) and X' — (a e ,a' h ) (same electron 
spin) and V^ a ' h is given by equation (|67l) . 

B. Qualitative description of the spin relaxation 
processes 

In order to describe qualitatively the variety of dif- 
ferent processes accounted for by Hamiltonian (1701) it is 
convenient to consider an initial state ip n as a linear com- 
bination of a dominant component \M)\ l) e \ f|")/j plus a 
minor contribution of two dark components, which arise 
from the coherent exchange of the Mn with either the 
electron or the hole: 

= \M)\i) e \it) h + 

+ e e \M - 1)| t) e | 1t)h + e h \M + 1)| |) e | ^72) 

where e e cx j e /jeh and €h oc jh/ jeh are small dimension- 
less coefficients that can be obtained doing perturbation 
theory. 

Depending on the elementary process that takes place, 
there are several possible final states: 

1. Hole spin relaxation. In this case the final state 
would be dominantly a dark exciton whose wave 
function read: 

|W) = \M)\l)e\V)h + 0(e) (73) 

and the scattering rate Tq would be proportional 
to |i^| 2 . This is process II in figure ©. 

2. Hole spin relaxation plus coherent hole-Mn spin 
flip. This is process III in figure (jB]) This can be re- 
alized through 2 dominant channels. An incoherent 
hole spin flip will couple the dominant component 
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FIG. 7: (Color online) Calculated rates for the transitions 
between exciton states in a Mn-doped quantum dot due to 
hole-phonon coupling. Left panel, red line, transition from 

— 1, — | > to | + 2, — I >, Right panel, red line, transition 
from| — 1, — | > to | — 2, — i >. blue line, transition from 

— 1, — | > to | + 1, — | >. L y is fixed at 6nm and scanning 
La, changes the LH-HH mixing parameter Jh±- 

of the initial state, \M)\ \.} e \ with a secondary 
component \M)\ l) e \ JJ.)/, of the final state 

= \M-l)\i) e \it) h + 
+ e h \M)\i) e \^) h + 0(e e ) (74) 

In this case the final state is a bright exciton in the 
same branch +1 than the initial state but the Mn 
component goes from M to M - 1. 

The second channel comes from the hole spin flip 
of the minority dark component of the initial state, 
€h\M + 1)| J,)e] ^}h which decays into the majority 
component of the final state 

|</v) = |Af + l)U) e |^) h + 0(e eifc ) (75) 

Thus, in this second case a hole spin flips due to 
phonons, plus a coherent Mn-electron spin flip con- 
nect the X = +1,M initial state to the X = 
— 1,M + 1 state. Thus, both the initial and final 
state in this process are the same than in the first 
channel, the rates for each would be proportional 
to e^To, but the decay pathways are different, and 
interferences are expected. 

3. Hole spin relaxation plus coherent electron-Mn spin 
flip. This is process I in figure © As in the previous 
case, there are two channels for this type of process. 
In the first channel, the majority component of the 
initial state decays into a final state given by: 

|W) = \M-l)\t)e\H)h + e' e \M)\l) c U)>i (76) 

The incoherent hole spin flip connects the initial 
state (fT2"|) to the final state (|7fJ)l through the mi- 
nority component |M)| \.) e \ jl)^ of the later. 
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FIG. 8: (Color online) Simulation of PL intensity from state 
(X — -\-l,M z — — |) under the influence of a driving laser 
pumping the system resonantly from optical ground state 
M K = -§ to the excited state (X = -1, M z = -§) for 
two laser intensities. The inset is the PL spectrum assum- 
ing all the states are equally populated. In the calculation, 
the quantum dot anisotropy (L;c=5nni and L a =6nm) controls 
the LH-HH mixing. The other parameters are discussed in the 
text. 

The second channel comes from the hole spin flip 
of the minority dark component of the initial state, 
e e \M — 1)| f) e | -fr)/j which decays into the majority 
component of the final state 

|W) = |M-l)|t>e|^>fc + 0(Ce,h) (77) 

Thus, a hole spin flips due the phonon, plus a coher- 
ent Mn-electron spin flip connect the X = +1,M 
initial state to the X = — 1, M — 1 state. The scat- 
tering rate of these two process scales as e^To 

C. Calculation of the relaxation rates 

In order to implement equations (|70|71|) to compute 
scattering rates, we use the single particle basis for the 
holes done with equations (|T2|) which leads, at finite mag- 
netic field to the matrix element that would be incor- 
porated into equations (J7T]) to compute the rates using 
equation (|B7I) . As discussed above, a zero field model 
(|12p yields a zero spin-flip matrix element in equation 
This is a feature of the simple hole model rather 
than an intrinsic property of the system. Thus, for the 
sake of simplicity, we compute the rates between exci- 
ton states by computing the matrix element (|68[) as if 
there was a magnetic field that yields the energy split- 
ting between the initial and final exciton states equal to 
the splitting produced by the exchange interaction with 
the Mn spin. 

In the calculation of the rates we perform an additional 
approximation: we only consider spin-flip terms in equa- 
tion (1711) and we do exclude spin-conserving terms. The 
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FIG. 9: (Color online) Left panel: average magnetization and 
right panel: occupation of the different spin states under op- 
tical pumping of the state (X = — 1, M z = — |). Parameters 
are the same as for the calculation presented in figure J5J. 



results for transition rates from the state n with domi- 
nant (—1,-1) to 3 possible final states with dominant 
components (+2, — ~), (— 2, — |) and (+1, — |) as a func- 
tion of the spin-flip Mn hole exchange Jh±_ , are shown in 
figure ([7]). The transition to the (+2,— |), which only 
involves the irreversible spin flip of the hole via a phonon 
emission is the dominant process and has a lifetime of 30 
ns. The transition to the (+1,— |) state requires both 
the hole spin flip and the Mn-hole spin flip and it is 3 
orders of magnitude less efficient. 

Thus, these calculations indicate that the most likely 
mechanism for Mn spin orientation in the presence of 
an exciton combines a rapid bright-to dark conversion, 
produced by phonon induced hole spin flip, and a dark to 
ground transition, enabled by Mn-carrier spin exchange 
and radiative recombination. 



VI. LASER DRIVEN SPIN DYNAMICS 

A. Summary of scattering mechanisms and master 
equation 

The spin dynamics of a single Mn atom in a laser driven 
quantum dot is described in terms of the 24 exciton states 
and the 6 ground states <p m . In the previous sections 
we have calculated the scattering rates of these states. 
They can be summarized as follows: 

1. Transitions from the ^ n to the <f> m , via photon 
emission (eq. (|41|l). In the case of bright excitons, 
this process is the quickest of all, with a typical life- 
time of 0.3 ns. In the case of dark excitons the life- 
time depends on the bright/dark mixing, which is 
both level and dot dependent. Dark lifetime ranges 
from twice the one of bright excitons to 1000 times 



larger, ie, between 1 and 300 nanoseconds. In any 
event, dark recombination involves a Mn spin flip. 

2. Transitions between different <fi m states, due to Mn 
spin phonon coupling (eq. (|53l) h The lifetimes of 
these transitions are, at least, 1ms (see right panel 
of figure [4|). 

3. Transitions between different exciton states 
that flip the spin of the Mn only, due to Mn-phonon 
coupling (eq. (|59|) . The lifetimes of these transi- 
tions are, at least, 0.1ms (see left panel figure QJ. 

4. Transitions between exciton states due to hole- 
phonon coupling (eq. (|70j) . The bright to dark 
transition is the quickest process with a lifetime of 
about 30ns (see figure [7]) . Bright to bright tran- 
sitions, combining holc-phonon and Mn-carrier in- 
teractions, have lifetimes in the 10 fj,s range. 

In addition to these dissipative scattering processes, 
we have to consider driving effect of the laser field, de- 
scribed in the semiclassical approximation. All things 
considered, we arrive to a master equation that describes 
the evolution of the occupations pn, where N = (n, m) 
includes states both with and without exciton in the dot. 
The master equation reads: 



0.015 



n^n'Pn 



(78) 
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N> 



Eq. (|78|) is a system of 36 coupled differential equations 
that we solve by numerical iteration, starting from a ther- 
mal distribution for the initial occupation p^. Since the 
temperature is larger than the energy splitting in the 
ground state, but much smaller than the band gap, at 
t = we have the six ground states with similar occu- 
pation P m ~ 1/6, P n ~ 0. As a result, the average 
magnetization, defined as: 

(M z ) = Y,Pm(<t>m\M z \4> m ) (79) 
m 

is zero, at zero magnetic field, as expected. 

B. Optical Mn spin orientation 

Under the action of the laser, the exciton states be- 
come populated and, under the adequate pumping con- 
ditions, the average Mn magnetization (M z ) acquires a 
non-zero value. This transfer of angular momentum, 
known as optical Mn spin orientation has been observed 
experimentally^ and predicted theoretically^. It results 
from a decrease of the Mn spin lifetime in the pres- 
ence of the exciton in the dot. In that circumstance, 
the laser transfer population from the M z state to the 
X, M z state. The enhanced relaxation transfers popula- 
tion from X, M z to X, M' z and the recombination to M' z 
state. Thus, if the laser is resonant with a single M z to 
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FIG. 10: (Color online) Evolution of Mn spin orientation 
efficiency as a function of the laser power. The pumping is 
detected in the PL intensity from state (X = +1,M Z = — |) 
under the influence of a driving laser pumping the system from 
optical ground state M z = — | to the excited state (X = — 1, 
Mz = — §)• The inset presents the laser power dependence 
of Tp um p, from where we can see that the efficiency of the 
pumping gets higher with the increasing of the laser power. 



Quantity 


Symbol 


Value 


Hole-Mn exchange 


h 


0.31 meV 


Electron-Mn exchange 




-0.09 meV 


Electron-Hole 


jeh 


-0.73 meV 


Uniaxial Anisotropy 


D 


10 fieV 


In plane Anisotropy 


E 





Quantum dot width 


Ly 


6nm 


Quantum dot width 




5nm 


Quantum dot height 




3nm 



TABLE I: Parameters used in the simulation of the resonant 
PL observed in the time resolved optical pumping experi- 
ments. 



X, M z transition, the M z state is depleted, which results 
in a decrease of the PL coming both from the X, M z and 
the —X, M z transitions. 

In figure (jHJ we show the result of our simulations for 
a dot at thermal equilibrium (fc^T = 4K) at t = 
which is pumped with a laser pulse resonant with the 
X = —1,M Z = — | transition, which is the high en- 
ergy one, since the hole is parallel to Mn spin. The 
laser pulse has a duration of 300 nanoseconds, so that 
the spectral broadening is negligible. In the upper panel 
we plot the PL coming from the counter-polarized transi- 
tion, X = +1, M z = — |, which has lower energy and can 
be detected without interference with the laser, for two 
different pumping power intensity. It is apparent that af- 
ter a rise of the PL in a time scale of 8 ns, corresponding 
the spin relaxation of the exciton spin from X = — 1 to 
X = +1, the PL signal is depleted. The origin of the 
depletion is seen in figure ©■ The occupation of the 
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Time(ns) 

FIG. 11: (Color online) Evolution of Mn spin orientation effi- 
ciency as a function of the valence band mixing controlled by 
the anisotropy of the confinement potentiel (L H = 13nm, vari- 
able La,). The pumping is detected in the PL intensity from 
state (X = +1,M% = — §) under the influence of a driving 
laser pumping the system from optical ground state M z = — | 
to the excited state (X = —1, M z = — |). The inset shows the 
evolution of T pump with L x and Jh± ■ The exciton generation 
rate is fixed at 1 ns _1 . 



M z = — | spin state in the ground reduced down to zero, 
in benefit of the other Mn spin states. 

Accordingly, the average magnetization becomes finite. 
Thus, net angular momentum is transferred from the 
laser to the Mn spin. The transfer takes place through 
Mn spin relaxation enabled in the presence of the exci- 
ton. As discussed above, the most efficient mechanism 
combines hole-spin relaxation due to phonons combined 
with dark-bright mixing, which involves a Mn spin flip. 

Interestingly, the fact that in the steady state several 
Mn spin states are occupied, including the higher energy 
ones, is compatible with a picture in which the Mn spin 
is precessing. Thus, a steady supply of spin-polarized 
excitons in the dot would result in the precession of the 
Mn spin, a scenario similar to that of current drive spin- 
torque oscillators 45 . Further work necessary to confirm 
this scenario is outside the scope of this paper. 

The efficiency of the process increases with the laser 
power, as shown in figure (fTOl) . We define the spin ori- 
entation time Tp Um p as the time at which the PL of the 
counter polarized transition is half the maximum. We 
can see that, as expected, Tp Ump is a decreasing func- 
tion of the laser power. A pumping time Tp Ump ~90ns 
is obtained with a generation rate of about Ins -1 . The 
amplitude of the valence band mixing, controlled by the 
anisotropy of the confinement potential or the in-plane 
strain distribution, is the main quantum dot parame- 
ter controlling the efficiency of the optical pumping. As 
presented in figure (|11[) . decreasing the quantum dot 
anisotropy, i.e., decreasing the LH-HH mixing parameter 



J/j_l, produces a rapid increase of r pump (inset of figure 
(fTT|) . This is a direct consequence of the reduction of the 
phonon induced hole spin flip. 



VII. SUMMARY AND CONCLUSIONS 

We have studied the spin dynamics of a single Mn atom 
in a CdTe quantum dot excited by a laser that drives the 
transition between the 6 optical ground states, associ- 
ated to the 25* + 1 states of the Mn spin S = | , and the 
24 single exciton states, corresponding to X = ±1,±2 
states interacting with the Mn spin. The main goal is 
to have a microscopic theory for the Mn spin relaxation 
mechanisms that makes it possible to produce laser in- 
duced Mn spin orientation in a time scale of less than 
100 ns»i _ — For that matter, we need to describe how the 
Mn and the quantum dot exciton affect each other. 

In section (pTj) we describe the different terms in the 
Mn spin Hamiltonian, including exchange with the 0- 
dimensional exciton. The symmetry of the exchange in- 
teraction depends on the spin properties of the carriers, 
which in the case of holes are strongly affected by the 
interplay of confinement, strain and spin orbit coupling. 
In section (JTTJ) we also use a model for hole a 18 ' 23 i 30 in 
quantum dots, which permits to obtain analytical ex- 
pressions for the wave functions of the holes, the hole-Mn 
exchange, in terms of the dimensions of the dot and the 
Kohh-Luttinger Hamiltonian. 

In section (IIIII) we study the dissipative dynamics of 
the Mn spin due to its coupling to phonons, both with 
and without excitons in the dot. The Mn spin-phonon 
coupling arises from the time dependent stochastic fluc- 
tuations of the crystal field and thereby of the single 
ion magnetic anisotropy, induced by the phonon field. 
Whereas the Mn spin relaxation is accelerated by 2 
or 3 orders of magnitude in the presence of the exci- 
ton, the efficiency of this mechanism is too low to ac- 
count for the optical orientation of the Mn spin reported 
experimentally-^—. The small Mn spin-phonon coupling 
comes from the small magnetic anisotropy of Mn as a 
substituional impurity in CdTe. 

In section pV[) we describe the interaction between 
the hole spin and the phonons in non-magnetic dots. Us- 
ing the simple analytical model for the holes presented 
in section |TT] we obtain analytical formulas for the hole 
spin relaxation. We find that hole spin lifetime can be 
in the range of 30 ns for a hole spin splitting as large 
as that provided by the hole-Mn coupling. Thus, we ex- 
pect that bright excitons will relax into dark excitons 
via hole-spin relaxation. This provides a microscopic 
mechanism to the scenario for Mn spin relaxation pro- 
posed by Cywinski^: bright excitons relax into dark ex- 
citons, via carrier spin relaxation, and the joint process 
of Mn-carrier spin exchange couples the dark excitons 
to the bright excitons, resulting in PL from dark states 
which implies Mn spin relaxation in a time scale of a few 
nanoseconds. This scenario is confirmed by calculations 
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presented in section fVl Finally, in section |VT1 we present 
the master equation that governs the dynamics of the 30 
states of the dot, we solve it numerically and we model 
the optical Mn spin orientation reported experimentally. 
Our main conclusions are: 

• Mn spin-phonon spin relaxation is presumably too 
weak to account for Mn spin dynamics in the pres- 
ence of the exciton 

• The Mn spin orientation is possible in a time scale 
of one hundred nanoseconds via a combination of 
phonon-induced hole spin relaxation and the sub- 
sequent recombination of the dark exciton enabled 
by spin-flip exchange of the Mn and the carrier 

• The critical property that governs the hole-Mn ex- 
change and the hole spin relaxation is the mixing 
between light and heavy holes, which depends both 
on the shape of the dot and on strain. 

• Our microscopic model permits to account for the 
optically induced Mn spin orientation. 

Future work should address how the coupling of the 
electronic Mn spin to the nuclear spin modifies our re- 



sults. This probably plays a role for the Mn in the dot 
without excitons. In addition, future work should study 
the role played by Mn spin coherence, and the interplay 
between optical and spin coherence. 
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Appendix A: Kohn Luttinger Hamiltonian 

The Kohn-Luttinger Hamiltonian for the 4 topmost 
valence bands of a Zinc Blend compound are given by: 



T~L(k Xl ky, k z ) — 



/ P + Q-1kv b B 
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and 



(A4) 



where 71,2,3 are dimensionless material dependent pa- 
rameters, 7 = |(72 + 73), /•* = 5(72 - 73), m is the free 
electron mass, k\ — k 2 . + k 2 and k± = k x ± ik y . For the 
dot states the relevant parameters are: 
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(A5) 



Appendix B: General formula for phonon-induced 
spin-flip rate 

In this appendix we derive a general formula for the 
scattering rate between two electronic state n and n' in- 
duced by a phonon emission. The Hamiltonian of the 
system can be split in 3 parts, the electronic states n, the 
phonon states, and their mutual coupling. The phonon 
states are labelled according to their polarization and 
momentum, A, q. We consider the following coupling 



V 



E 

m,m' ,q,\ 



\m)(m'\{b{ q 



bx.-. 



(Bl) 
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where m and m! are electronic states. We refer to the 
free phonon states as the reservoir states. Within the 



J 



Born-Markov approximation, the scattering rate between 
states n and n' . 



IW =|E P ^E \(nr\V\n'r')f S (E n - E n , + e r - e r ,) 



(B2) 



where P r is the occupation of the r reservoir state with 
energy e r . This equation can be interpreted as a statis- 
tical average over reservoir initial states r of the Fermi 
Golden rule decay rate of state N, r. 

The sums over r and r' are performed using the fol- 
lowing trick. For a given r, the initial reservoir state, r' 
must have an additional phonon, since we consider the 
phonon emission case. Thus, we write: 



10 



1 



>/n\', q > + 1 



so that 



(B3) 



Wlx + b- q ,x\r') = <U^a,v + 1 (B4) 



The matrix element 



(nr\V\nV) = V? A " Jn\,, q , + 1 



(B5) 



We see how from all the terms in the sum that defines the 
coupling, only one survives and fixes the index r'. Thus, 
the only the sums left are the over the initial reservoir 
states and the A, q index that define the final state. Now 
we use the definition of the Bose function: 

J2 Pr K, 9 < + 1) = n B (w A , (</)) + 1 (B6) 

r 

and we arrive to the following expression for the rate: 



(B7) 



Notice that it is possible to write the rate as a sum 
over different contributions arising from different po- 
larizations, r = X^a-^- I n t ne particular case that 
we can neglect the dependence of the matrix element 
Vq,\(n,n') ~ |V( n! „')| 2 on q and A, we arrive to the fol- 
lowing expression: 



2-7T 

= — (n B (A) + l)p A (A)|V n>r , 



(B8) 



where p(A) = ^ A 6 (A — tj\(q)) is the density of states 
of the phonons evaluated at the transition energy A. 
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